Strong exciton-plasmon coupling in semiconducting carbon nanotubes 
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We study theoretically the interactions of excitonic states with surface electromagnetic modes 
of small-diameter ( < 1 nm) semiconducting single-walled carbon nanotubes. We show that these 
interactions can result in strong exciton-surface-plasmon coupling. The exciton absorption line 
shape exhibits Rabi splitting ~ 0.1 eV as the exciton energy is tuned to the nearest interband surface 
plasmon resonance of the nanotube. We also show that the quantum confined Stark effect may be 
used as a tool to control the exciton binding energy and the nanotube band gap in carbon nanotubes 
in order, e. g., to bring the exciton total energy in resonance with the nearest interband plasmon 
mode. The exciton-plasmon Rabi splitting we predict here for an individual carbon nanotube is 
close in its magnitude to that previously reported for hybrid plasmonic nanostructures artificially 
fabricated of organic semiconductors on metallic films. We expect this effect to open up paths to 
new tunable optoelectronic device applications of semiconducting carbon nanotubes. 



PACS numbers: 78.40.Ri, 73.22.-f, 73.63.Fg, 78.67.Ch 



I. INTRODUCTION 

Single-walled carbon nanotubes (CNs) are quasi-one- 
dimensional (ID) cylindrical wires consisting of graphene 
sheets rolled- up into cylinders with diameters ~ 1 — 10 nm 
and lengths - 1 - 10 4 fim 0, S, S, CNs are 

shown to be useful as miniaturized electronic, electrome- 
chanical, and chemical devices [f|, scanning probe de- 
vices [6j, and nanomaterials for macroscopic compos- 
ites 0] ■ The area of their potential applications was re- 
cently expanded to nanophotonics d, 0, EH E3 after the 
demonstration of controllable single-atom incapsulation 
into CNs [HI, 0, HH , and even to quantum cryptog- 
raphy since the experimental evidence was reported for 
quantum correlations in the photoluminescence spectra 
of individual nanotubes [l6| . 

For pristine (undoped) single-walled CNs, the numer- 
ical calculations predicting large exciton binding ener- 
gies (~ 0.3 - 0.6 eV) in semiconducting CNs [H, 
19( | and even in some small-diameter (~ 0.5 nm) metal- 
lic CNs [2(J, followed by the results of various exciton 
photoluminescence measurements [HI, [U, [52, [H, [24| [25[ , 
have become available. These works, together with other 
reports investigating the role of effects such as intrinsic 
defects [H, [26| , exciton-phonon interactions [24l [261 . l27l . 
Hi Hil , external magnetic and electric fields |30l. l3ll. l32j] . 
reveal the variety and complexity of the intrinsic optical 
properties of CNs [33l ] . 

Here we develop a theory for the interactions between 
excitonic states and surface electromagnetic (EM) modes 
in small-diameter (< 1 nm) semiconducting single- walled 
CNs. We demonstrate that such interactions can result 
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in a strong exciton-surface-plasmon coupling due to the 
presence of low-energy (~ 0.5 — 2 eV) weakly-dispersive 
interband plasmon mod es [34l| and large exciton excita- 
tion energies ~ 1 eV (35l.l36l| in small-diameter CNs. Pre- 
vious studies have been focused on artificially fabricated 
hybrid plasmonic nanostructures, such as dye molecules 
in organic polymers deposited on metallic films [37| , semi- 
conductor quantum dots coupled to metallic nanopar- 
ticles 38 1, or nanowires [39l |. where one material car- 
ries the exciton and another one carries the plasmon. 
Our results are particularly interesting since they reveal 
the fundamental EM phenomenon — the strong exciton- 
plasmon coupling — in an individual quasi-lD nanos- 
tructure, a carbon nanotube. 

The paper is organized as follows. Section [IT] presents 
the general Hamiltonian of the exciton interaction with 
vacuum-type quantized surface EM modes of a single- 
walled CN. No external EM field is assumed to be ap- 
plied. The vacuum-type-field we consider is created by 
CN surface EM fluctuations. In describing the exciton- 
field interaction on the CN surface, we use our re- 
cently developed Green function formalism to quantize 
the EM field in the presence of quasi- ID absorbing bod- 
ies 0, SI E2, S H El- The formalism follows the 
original line of the macroscopic quantum electrodynamics 
(QED) approach developed by Welsch and coworkers to 
correctly describe medium-assisted electroma gne tic vac- 
uum effects in dispersing and absorbing media [4rj l47l lis]] 
(also refs. therein). Section HTT1 explains how the inter- 
action introduced in Sec. II results in the coupling of 
the excitonic states to the nanotube's surface plasmon 
modes. Here we derive, calculate and discuss the charac- 
teristics of the coupled exciton-plasmon excitations, such 
as the dispersion relation, the plasmon density of states 
(DOS), and the optical absorption line shape, for par- 
ticular semiconducting CNs of different diameters. We 
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quantized field Hamiltonian is 



FIG. 1: (Color online) The geometry of the problem. 



also analyze how the electrostatic field applied perpen- 
dicular to the CN axis affects the CN band gap, the ex- 
citon binding energy, and the surface plasmon energy, 
to explore the tunability of the exciton-surface-plasmon 
coupling in CNs. The summary and conclusions of the 
work are given in Sec. lIVI All the technical details about 
the construction and diagonalization of the exciton-field 
Hamiltonian, the EM field Green tensor derivation, the 
perpendicular electrostatic field effect, are presented in 
the Appendices in order not to interrupt the flow of the 
arguments and results. 



II. EXCITON-ELECTROMAGNETIC-FIELD 
INTERACTION ON THE NANOTUBE SURFACE 

We consider the vacuum-type EM interaction of an ex- 
citon with the quantized surface electromagnetic fluctua- 
tions of a single- walled semiconducting CN by using our 
recently developed Green function formalism to quan- 
tize the EM field in thepresence of quasi- ID absorbing 
bodies [H El, E3, H| 0, |||. No external EM field 
is assumed to be applied. The nanotube is modelled by 
an infinitely thin, infinitely long, anisotropically conduct- 
ing cylinder with its surface conductivity obtained from 
the realistic band structure of a particular CN. Since the 
problem has the cylindrical symmetry, the orthonormal 
cylindrical basis {e r , e^, e 2 } is used with the vector e z di- 
rected along the nanotube axis as shown in Fig. [T] Only 
the axial conductivity, a zz , is taken into account, whereas 
the azimuthal one, a vv , being strongly suppressed by the 
transverse depolarization effect [49l.l50l.l5ll. l52l. l53l. l54j . is 
neglected. 

The total Hamiltonian of the coupled exciton-photon 
system on the nanotube surface is of the form 



H = Hi 



He 



Hi, 



(1) 



where the three terms represent the free (medium- 
assisted) EM field, the free (non-interacting) exciton, and 
their interaction, respectively. More explicitly, the second 
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where the scalar bosonic field operators /^(n,w) and 
/(n, uj) create and annihilate, respectively, the surface 
EM excitation of frequency lj at an arbitrary point 
n = R„ = {Rcn, tfn, z n } associated with a carbon atom 
(representing a lattice site - Fig. [I} on the surface of the 
CN of radius Rcn- The summation is made over all the 
carbon atoms, and in the following it is replaced by the 
integration over the entire nanotube surface according to 
the rule 
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where S = (3\/3/4)& 2 is the area of an elementary equi- 
lateral triangle selected around each carbon atom in a 
way to cover the entire surface of the nanotube, b= 1.42 A 
is the carbon-carbon interatomic distance. 

The second quantized Hamiltonian of the free exciton 
(see, e.g., Ref. (55[) on the CN surface is of the form 

H ex = J2Ef(n)Bt +mJ B mJ = Y,E f (k)Bl f B k j, (4) 

n,m,/ k,/ 

where the operators B^ i and B n j create and annihi- 
late, respectively, an exciton with the energy Ef(n) in 
the lattice site n of the CN surface. The index / 0) 
refers to the internal degrees of freedom of the exciton. 
Alternatively, 

B L = ^E</ e4kn and ^,/ = K / ) t (5) 



create and annihilate the /-internal-state exciton with 
the quasi- momentum k = {fc l/ ,,/c z }, where the azimuthal 
component is quantized due to the transverse confine- 
ment effect and the longitudinal one is continuous, N is 
the total number of the lattice sites (carbon atoms) on 
the CN surface. The exciton total energy is then written 
in the form 



E f (k)=EiH(k, 



2 Ul 



h k 



(6) 



v ' 2M ex {k v ) 
Here, the first term represents the excitation energy 

EiH(k v )=E g (k v ) + Ej ) f \k v ) (7) 

of the /-internal-state exciton with the (negative) bind- 
( f) 

ing energy E^ , created via the intcrband transition with 
the band gap 



E g {k {p ) = e e (k v ) + e h (kp), 



(8) 



where e et h are transversely quantized azimuthal electron- 
hole subbands (see the schematic in Fig. \2§ ■ The second 
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FIG. 2: (Color online) Schematic of the two transversely 
quantized azimuthal electron-hole subbands (left), and the 
first-interband ground- internal-state exciton energy (right) in 
a small-diameter semiconducting carbon nanotube. Subbands 
with indices j — 1 and 2 are shown, along with the optically 
allowed (exciton-related) interband transitions See text 
for notations. 



being the interaction matrix element where the exciton 
with the energy E ex ' c = Tuvf is excited through the elec- 
tric dipole transition (d£) z = (0|(d n ) a |/) in the lattice 
site n by the nanotube's transversely (longitudinally) po- 
larized surface EM modes. The modes are represented in 
the matrix element by the transverse (longitudinal) part 
of the Green tensor zz-component G zz (n, m,oj) of the 
EM subsystem (Appendix B). This is the only Green ten- 
sor component we have to take into account. All the other 
components can be safely neglected as they are greatly 
suppressed by the strong transverse depolarization effect 
in CNs [H, 0, HH [H, HI, As a consequence, only 
<Jzz(Rcn the axial dynamic surface conductivity per 
unit length, is present in Ea.(|12p. 

Equations ([I])-(|12p form the complete set of equations 
describing the exciton-photon coupled system on the CN 
surface in terms of the EM field Green tensor and the 
CN surface axial conductivity. 

III. EXCITON-SURFACE-PLASMON 
COUPLING 



term in Eq. ^ represents the kinetic energy of the trans- 
lational longitudinal movement of the exciton with the ef- 
fective mass M ex = m e + mn, where m e and are the 
(subband-dependent) electron and hole effective masses, 
respectively. The two equivalent free-exciton Hamilto- 
nian representations are related to one another via the 
obvious orthogonality relationships 

1 £ e ^( k - k ').n =5kk() ^£ e - i( n- m ). k = (Wn (g) 
n k 

with the k-summation running over the first Brillouin 
zone of the nanotube. The bosonic field operators in Hp 
are transformed to the k-representation in the same way. 

The most general (non-relativistic, electric dipole) 
exciton-photon interaction on the nanotube surface can 
be written in the form (we use the Gaussian system of 
units and the Coulomb gauge; see details in Appendix IA")) 

/>oo 

H int =J2 / <Mg/" ) (n,m,w)B£ / 

n,mJ J ° 

- gjr 5 (n, m, w)B nJ ] /(m, u) + h.c, (10) 

where 

g ( , ±) (n,m,w) = g|(n,m,w)± — g'i(n,m,w) (11) 
with 



g| (ll) (n,m,w) = -i—^-^/nfiuj Re a zz (R CN , u) 

x (d f n ) z u ^G zz (n,m,uj) (12) 



For the following it is important to realize that the 
transversely polarized surface EM mode contribution to 
the interaction (fTT))) - (fT2")) is negligible compared to the 
longitudinally polarized surface EM mode contribution. 
As a matter of fact, ± G ZZ (n, m, u>) =0 in the model of an 
infinitely thin cylinder we use here (Appendix [B]) , thus 
yielding 

gl(n,m,w) = ! g ( f ±) (n,m,w)=±— g''(n,m,tj) (13) 

in Eqs. (fTtI)) - (fT2"|) . The point is that, because of the 
nanotube quasi-one-dimensionality, the exciton quasi- 
momentum vector and all the relevant vectorial matrix 
elements of the momentum and dipole moment operators 
are directed predominantly along the CN axis (the lon- 
gitudinal exciton; see, however, Ref. [56|). This prevents 
the exciton from the electric dipole coupling to trans- 
versely polarized surface EM modes as they propagate 
predominantly along the CN axis with their electric vec- 
tors orthogonal to the propagation direction. The longi- 
tudinally polarized surface EM modes are generated by 
the electronic Coulomb potential (see, e.g., Ref. [57|). 
and therefore represent the CN surface plasmon excita- 
tions. These have their electric vectors directed along the 
propagation direction. They do couple to the longitudi- 
nal excitons on the CN surface. Such modes were ob- 
served in Ref. [HJ|. They occur in CNs both at high 
energies (well-known 7r-plasmon at ~ 6 eV) and at com- 
paratively low energies of ~ 0.5 — 2 eV. The latter ones 
are related to the transversely quantized interband (inter- 
van Hove) electronic transitions. These weakly-dispersive 
modes [13, HI| are similar to the intersubband plasmons 
in quantum wells [5^ |. They occur in the same energy 
range of ~ 1 eV where the exciton excitation energies 
are located in small-diameter (< 1 nm) semiconducting 
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CNs [35|, [36J . In what follows we focus our consideration 
on the exciton interactions with these particular surface 
plasmon modes. 



A. The dispersion relation 

To obtain the dispersion relation of the coupled exci- 
ton- plasmon excitations, we transfer the total Hamil- 
tonian (|TJ)— (JTOj) and (fl~3"|) to the k-representation using 
Eqs. (O and ([5]), and then diagonalize it exactly by means 
of Bogoliubov's canonical transformation technique (see, 
e.g., Ref. The details of the procedure are given in 

Appendix [Cj The Hamiltonian takes the form 

H = £ ^(k) (1 (k)4(k) + E . (14) 

k,/j=l,2 

Here, the new operator 

4(k) = [<(k,w/)Bk,/ -^(k,w/)Bi ki/ ] (15) 
/ 

+ Jlu [^(k, w)/(k, u) - v;(k, «)/t(-k, La) 

annihilates and £t(k)== [^(k)]* creates the exciton-plas- 
mon excitation of branch /z, the quantities and 
are appropriately chosen canonical transformation coef- 
ficients. The "vacuum" energy E represents the state 
with no exciton-plasmons excited in the system, and 
hcu^h) is the exciton-plasmon energy given by the solu- 
tion of the following (dimcnsionless) dispersion relation 



. °° xTj{x)p{x) 
'-f- £ f - / dx ~ ~2 ZT- 

n .10 X u ~~ X 



(16) 
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FIG. 3: (Color online) (a),(b) Calculated dimensionless (see 
text) axial surface conductivities for the (11,0) and (10,0) 
CNs. The dimensionless energy is defined as [Energy] /2 r yo, 
according to Eq. (|17[) . 
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with 7o = 2.7 eV being the carbon nearest neighbor over- 
lap integral entering the CN surface axial conductivity 
<Jzz{Rcn, <jj). The function 



4\4\ 2 x 3 
3hc 3 



2to 

h 



(18) 



with d{ — 5Z n (0|(d n )z|/) represents the (dimensionless) 
spontaneous decay rate, and 



p(x) 



3So 



16iraR CN 



Re- 



1 



(19) 



stands for the surface plasmon density of states (DOS) 
which is responsible for the exciton decay rate varia- 
tion due to its coupling to the plasmon modes. Here, 
a = e 2 /hc = 1/137 is the fine-structure constant and 



a zz = 2irhcr zz /e 2 is the dimensionless CN surface axial 
conductivity per unit length. 

Note that the conductivity factor in Eq. (| 19[) equals 



Re 
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o zz [x] 



Im- 



1 



4ac 



Rcn \2^ x J e zz (x) - 1 



(20) 



in view of Eq. (|17[) and equation a zz — — iu)(e zz — l)/4irSpT 
representing the Drude relation for CNs, where e zz is the 
longitudinal (along the CN axis) dielectric function, S 
and pt are the surface area of the tubule and the number 
of tubules per unit volume, respectively [4lL 53. [50(| ■ This 
relates very closely the surface plasmon DOS func- 
tion (fl9|) to the loss function — Im(l/e) measured in Elec- 
tron Energy Loss Spectroscopy (EELS) experiments to 
determine the properties of collective electronic excita- 
tions in solids [341 ]. 

Figure [3] shows the low-energy behaviors of the func- 
tions a zz (x) and Re[l/CT zz (x)] for the (11,0) and (10,0) 
CNs (Rcn = 0.43 nm and 0.39 nm, respectively) we 
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FIG. 4: (Color online) (a),(b) Surface plasmon DOS and conductivities (left panels), and lowest bright exciton dispersion when 
coupled to plasmons (right panels) in (11,0) and (10,0) CNs, respectively. The dimensionless energy is defined as [Energy] /2 r yo, 
according to Eq. fTTJ). See text for the dimensionless quasi- momentum. 



study here. We obtained them numerically as follows. 
First, we adapt the nearest-neighbor non-orthogonal 
tight-binding approach [6l[ to determine the realistic 
band structure of each CN. Then, the room-temperature 
longitudinal dielectric functions e zz are calculated within 
the random-phase approximation [H,[63|], which are then 
converted into the conductivities a zz by means of the 
Drude relation. Electronic dissipation processes are in- 
cluded in our calculations within the relaxation-time ap- 
proximation (electron scattering length of 130Rcn was 
used We did not include excitonic many-electron 

correlations, however, as they mostly affect the real con- 
ductivity Re(<7 zz ) which is responsible for the CN optical 
absorption [l|| [2^, HU , whereas we are interested here 
in Re(l/(7 ZZ ) representing the surface plasmon DOS ac- 
cording to Eq. I|19p . This function is only non-zero when 
the two conditions, lm]a zz (x)] = and Re[a zz (cc)] — > 0, 
are fulfilled simultaneously [H, H^, [62j . These result in 
the peak structure of the function Re(l/cr zz ) as is seen in 



Fig. [3J It is also seen from the comparison of Fig. [3J (b) 
with Fig. [3J (a) that the peaks broaden as the CN di- 
ameter decreases. This is consistent with the stronger 
hybridization effects in smaller-diameter CNs [6J|. 

Left panels in Figs. H^a) and 0Jb) show the lowest- 
energy plasmon DOS resonances calculated for the (11,0) 
and (10,0) CNs as given by the function p{x) in Eq. (|19[) . 
Also shown there are the corresponding fragments of 
the functions Re[cr zz (:r)] and Im[tr zz (a:)]. In all graphs 
the lower dimensionless energy limits are set up to be 
equal to the lowest bright exciton excitation energy 

[E^c = 1.21 eV (x = 0.224) and 1.00 eV (x = 0.185) 
for the (11,0) and (10,0) CN, respectively, as reported in 
Ref.JIH by directly solving the Bethe-Salpeter equation]. 
Peaks in p(x) are seen to coincide in energy with zeros 
of Im[(T zz (a:)] {or zeros of Re[e zz (x)]}, clearly indicating 
the plasmonic nature of the CN surface excitations under 
consideration [58|,[65|]. They describe the surface plasmon 
modes associated with the transversely quantized inter- 
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band electronic transitions in CNs [58j]. As is seen in 
Fig. 0] (and in Fig. [3]), the interband plasmon excitations 
occur in CNs slightly above the first bright exciton ex- 
citation energy [53j, in the frequency domain where the 
imaginary conductivity (or the real dielectric function) 
changes its sign. This is a unique feature of the com- 
plex dielectric response function, the consequence of the 
general Kramers-Kronig relation (4(| . 

We further take advantage of the sharp peak structure 
of p(x) and solve the dispersion equation ([To]) for x^ an- 
alytically using the Lorentzian approximation 



p(x) 



p(x p )Axl 



(x — Xp) 2 + Ax 2 



(21) 



Here, x v and Ax p are, respectively, the position and the 
half-width-at-half-maximum of the plasmon resonance 
closest to the lowest bright exciton excitation energy in 
the same nanotube (as shown in the left panels of Fig. [4} . 
The integral in Eq. ([TBI) then simplifies to the form 



fix 



xT f (x)p(x) F(x p )Ax 2 p 



Xp\ — X 2 



F(x p )Aa 



dx 



x fi — x p Jo ( x ~ x p) 2 + ^ x p 



arctan 



( x p \ + 1 
\Ax p J 2 



with F(x p ) = 2x p r{ ) (xp)p(xp)/'K. This expression is valid 
for all Xp apart from those located in the narrow inter- 
val (x p — Ax p , x p + Ax p ) in the vicinity of the plasmon 
resonance, provided that the resonance is sharp enough. 
Then, the dispersion equation becomes the biquadratic 
equation for Xp with the following two positive solutions 
(the dispersion curves) of interest to us 



Xl,2 



± 



Fps f 



(22) 



Here, F p = AF(xp)Axp{~K — Ax p /x p ) with the arctan- 
function expanded to linear terms in Ax p /x p <C 1. 

The dispersion curves (|22[) are shown in the right pan- 
els in Figs. |D(a) and 2Kb) as functions of the dimension- 
less longitudinal quasi-momentum. In these calculations, 
we estimated the interband transition matrix element in 
T f (x p ) [Eq. (JTHJ)] from the equation Id/1 2 = 3ftA 3 /4T™ d 
according to Hanamura's general theory of the exci- 
ton radiative decay in spatially confined systems [6t| . 
where r™ d is the exciton intrinsic radiative lifetime, and 
A = 2irch/E with E being the exciton total energy given 
in our case by Eq. ©. For zigzag-type CNs considered 
here, the first Brillouin zone of the longitudinal quasi- 
momentum is given by —2n~h/3b < fik z < 2irh/?>b [HQ. 
The total energy of the ground-intcrnal-state exciton can 
then be written as E = E exc + (2irh/3b) 2 t 2 /2M ex with 
— 1 < t < 1 representing the dimensionless longitudi- 
nal quasi-momentum. In our calculations we used the 
lowest bright exciton parameters E, 



(ii) 



1.21 eV and 



1.00 eV, T™ d = 14.3 ps and 19.1 ps, M ex = 0.44m and 
0.19mo (mo is the free-electron mass) for the (11.0) CN 
and (10,0) CN, respectively, as reported in Ref.[3a| by 
directly solving the Bethe-Salpeter equation. 

Both graphs in the right panels in Fig. 0] are seen to 
demonstrate a clear anticrossing behavior with the (Rabi) 
energy splitting ~ 0.1 eV. This indicates the formation 
of the strongly coupled surface plasmon-exciton excita- 
tions in the nanotubes under consideration. It is impor- 
tant to realize that here we deal with the strong exciton- 
plasmon interaction supported by an individual quasi- 
1D nanostructure — a single- walled (small-diameter) 
semiconducting carbon nanotube, as opposed to the ar- 
tificially fabricated metal-semiconductor nanostructures 
studied previuosly [3?], [H, H^] where the metallic com- 
ponent normally carries the plasmon and the semicon- 
ducting one carries the exciton. It is also important that 
the effect comes not only from the height but also from 
the width of the plasmon resonance as it is seen from the 
definition of the F p factor in Eq. (f2"2"| . In other words, 
as long as the plasmon resonance is sharp enough (which 
is always the case for interband plasmons), so that the 
Lorentzian approximation (I2ip applies, the effect is de- 
termined by the area under the plasmon peak in the DOS 
function (fl^)) rather than by the peak height as one would 
expect. 

However, the formation of the strongly coupled exci- 
ton-plasmon states is only possible if the exciton total en- 
ergy is in resonance with the energy of a surface plasmon 
mode. The exciton energy can be tuned to the nearest 
plasmon resonance in ways used for excitons in semicon- 
ductor quantum microcavities — thermally [o7l . l68l fsoj ] 
(by elevating sample temperature), or/and electrostati- 
cally [zO, El Iz2, uM ( y i a the quantum confined Stark 
effect with an external electrostatic field applied perpen- 
dicular to the CN axis). As is seen from Eqs. ([6]) and 
([7]) , the two possibilities influence the different degrees of 
freedom of the quasi- ID exciton — the (longitudinal) ki- 
netic energy and the excitation energy, respectively. Be- 
low we study the (less trivial) electrostatic field effect on 
the exciton excitation energy in carbon nanotubes. 



B. The perpendicular electrostatic field effect 

The optical properties of semiconducting CNs in an 
external electrostatic field directed along the nanotube 
axis were studied theoretically in Ref. [31(. Strong oscil- 
lations in the band-to-band absorption and the quadratic 
Stark shift of the exciton absorption peaks with the field 
increase, as well as the strong field dependence of the ex- 
citon ionization rate, were predicted for CNs of different 
diameters and chiralities. Here, we focus on the perpen- 
dicular electrostatic field orientation. We study how the 
electrostatic field applied perpendicular to the CN axis 
affects the CN band gap, the exciton binding/excitation 
energy, and the interband surface plasmon energy, to ex- 
plore the tunability of the strong exciton-plasmon cou- 
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pling effect predicted above. The problem is similar to the 
well-known quantum confined Stark effect first observed 
for the excitons in semiconductor quantum wells [7(J . 
However, the cylindrical surface symmetry of the exci- 
tonic states brings new peculiarities to the quantum con- 
fined Stark effect in CNs. In what follows we will gen- 
erally be interested only in the lowest internal energy 
(ground) excitonic state, and so the internal state index 
/ in Eqs. ((6|) and ||7J will be omitted for brevity. 

Because the nanotube is modelled by a continuous, in- 
finitely thin, anisotropically conducting cylinder in our 
macroscopic QED approach, the actual local symmetry 
of the excitonic wave function resulted from the graphene 
Brillouin zone structure is disregarded in our model (see, 
e.g., reviews (33l. l53jV The local symmetry is implicitly 
present in the surface axial conductivity though, which 
we calculate beforehand as described above [74| . 

We start with the Schrodinger equation for the elec- 
tron located at r e = {Rcn, fe, z e } and the hole located 
at Th = {RcNt iphi %h} on the nanotube surface. They 
interact with each other through the Coulomb potential 
V(r e ,r h ) = -e 2 /e|r e - r h \, where e = e zz (Q). The exter- 
nal electrostatic field F = {F, 0,0} is directed perpen- 
dicular to the CN axis (along the x-axis in Fig. [1} . The 
Schrodinger equation is of the form 



H e {F) + H h (F) + V(v e , v h )\ *(r e , r h ) = £*(r e , r h ) 

(23) 



2m h R 2 CN dif 2 h 



H 2 d 2 



2fi dz 2 



4>ex(z) = Eb4>ex(z), 



--Ehlphifh), 

(27) 
(28) 



where fi = m e mh/M ex is the exciton reduced mass, and 
Vgff is the effective longitudinal electron-hole Coulomb 
interaction potential given by 



2 r 2-rr r 2n 



V es( z ) =— / d(p e d<p h \ip e (ip e )\ 1^(^)1 V{tp e ,tp h ,z) 
e Jo Jo 

(29) 

with V being the original electron-hole Coulomb poten- 
tial written in the cylindrical coordinates as 



V(<p e ,<Ph,z) 



1 



{^+4^ iv sin 2 [(^-^)/2]}i/2- 



(30) 



with 



The exciton problem is now reduced to the ID equation 
(|28[) . where the exciton binding energy does depend on 
the perpendicular electrostatic field through the electron 
and hole subband functions ip e ,h given by the solutions of 
Eqs. (|2"6")l and (f2"T)) and entering the effective electron-hole 
Coulomb interaction potential ([29)) . 

The set of Eqs. (|2"5] ) -(|3"0" ]l is analyzed in Appendix O 
One of the main results obtained in there is that the 
effective Coulomb potential (f2"9")l can be approximated 
by an attractive cusp-type cutoff potential of the form 



d 2 



d 2 



2m e , h \R 2 CN dif 2 e h dz 2 e h 

(24) 

We further separate out the translational and relative 
degrees of freedom of the electron-hole pair by trans- 
forming the longitudinal (along the CN axis) motion 
of the pair into its center-of-mass coordinates given by 
Z = (m e z e + rrihZh) I 'M ex and z — z e — Zh- The exciton 
wave function is approximated as follows 



#(r e ,r h ) 



(l>ex(z)lp e ((pe)ll)h{tPh)- 



(25) 



The complex exponential describes the exciton center-of- 
mass motion with the longitudinal quasi-momentum k z 
along the CN axis. The function (f> ex (z) represents the 
longitudinal relative motion of the electron and the hole 
inside the exciton. The functions ip e ((p e ) and iph(fh) are 
the electron and hole subband wave functions, respec- 
tively, which represent their confined motion along the 
circumference of the cylindrical nanotube surface. 

Each of the functions is assumed to be normalized to 
unity. Equations (|23p and (f!M)) are then rewritten in view 
of Eqs. ©-© to yield 



d 2 



2m e R 2 CN dipl 



cRcnF cos (tp e 



(26) 



VeffW 



e[\z\+z (j,F)Y 



(31) 



where the cutoff parameter zq depends on the perpen- 
dicular electrostatic field strength and on the electron- 
hole azimuthal transverse quantization index j — 1, 2, ... 
(excitons are created in interband transitions involving 
valence and conduction subbands with the same quanti- 
zation index [53| as shown in Fig. [2]). Specifically, 



z (j,F)^2R CN 



ir - 21n2 [1 



A,(F)] 



7T + 21n2 [1 



A,(F)] 



(32) 



with Aj(F) given to the second order approximation in 
the electric field by 



Aj (F) w 2/j.M, 



ex 4 " j 

n 



(33) 



J l-2j l + 2j' 

where 9(x) is the unit step function. Approximation (f3Tj) 
is formally valid when zq (j, F) is much less than the exci- 
ton Bohr radius a* B (= eh 2 / fie 2 ) which is estimated to be 
~ IORcn for the first (j = 1 in our notations here) exciton 
in CNs (13. As is seen from Eqs. (|32|) and (f33|) . this is 
always the case for the first exciton for those fields where 
the perturbation theory applies, i. e. when Ai(F) < 1 in 
Eq. 
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Equation (|28j) with the potential (|3Tj) formally coin- 
cides with the one studied by Ogawa and Takagahara in 
their treatments of excitonic effects in fD semiconduc- 
tors with no external electrostatic field applied [76| . The 
only difference in our case is that our cutoff parameter 
(|32|) is field dependent. We therefore follow Ref. [76| and 

find the ground-state binding energy E^ 11 ^ for the first 
exciton we are interested in here from the transcendental 
equation 



In 



2z (l,F) 




0. (34) 



In doing so, we first find the exciton Rydberg energy, Ry* 
(= /j,e 4 /2h 2 e 2 ), from this equation at F — 0. We use the 
diameter- and chirality-dependent electron and hole ef- 
fective masses from Ref. [77j , and the first bright exciton 
binding energy of 0.76 eV for both (11,0) and (10,0) CN 
as reported in Ref. from ab initio calculations. We 
obtain Ry* = 4.02 eV and 0.57 eV for the (11,0) tube and 
(10,0) tube, respectively. The difference of about one or- 
der of magnitude reflects the fact that these are the semi- 
conducting CNs of different types — type-I and type-II, 
respectively, based on (2n + to) families [73}. The pa- 
rameters Ry* thus obtained are then used to find l-E^ 11 ^ 
as functions of F by numerically solving Eq. (f3"4"|) with 
z (l,F) given by Eqs. J35J) and ([33). 

The calculated (negative) binding energies are shown 
by the solid lines in Fig.[5Ja). Also shown there by dashed 
lines are the functions 



E^>{F) 



Ai(F)] 



(35) 



with Ai(F) given by Eq. (|33[) . They are seen to be fairly 
good analytical (quadratic in field) approximations to the 
numerical solutions of Eq. (|34|) in the range of not too 
large fields. The exciton binding energy decreases very 
rapidly in its absolute value as the field increases. Fields 
of only ~ 0.1 — 0.2 V//xm are required to decrease l-E^^I 
by a factor of ^2 for the CNs considered here. The rea- 
son is the perpendicular field shifts up the "bottom" of 
the effective potential (|3ip as shown in Fig. G^b) f° r the 
(11,0) CN. This makes the potential shallower and pushes 
bound excitonic levels up, thereby decreasing the exciton 
binding energy in its absolute value. As this takes place, 
the shape of the potential does not change, and the lon- 
gitudinal relative electron-hole motion remains finite at 
all times. As a consequence, no tunnel exciton ionization 
occurs in the perpendicular field, as opposed to the longi- 
tudinal electrostatic field (Franz-Keldysh) effect studied 
in Ref. [3l| where the non-zero field creates the potential 
barrier separating out the regions of finite and infinite 
relative motion and the exciton becomes ionized as the 
electron tunnels to infinity. 

The binding energy is only the part of the exciton ex- 
citation energy 0. Another part comes from the band 
gap energy ([5]), where e e and eh are given by the solu- 
tions of Eqs. (f2"6"]l and (j2"T|) , respectively. Solving them to 
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FIG. 5: (Color online) (a) Calculated binding energies of the 
first bright exciton in the (11,0) and (10,0) CNs as functions of 
the perpendicular electrostatic field applied. Solid lines are the 
numerical solutions to Eq. (I34|l , dashed lines are the quadratic 
approximations as given by Eq. (|35[) . (b) Field dependence of 
the effective cutoff Coulomb potential (|31|) in the (11,0) CN. 
The dimensionless energy is defined as [Energy]/ '270, accord- 
ing to Eq. (fT7)l . 



the leading (second) order perturbation theory approxi- 
mation in the field (Appendix iDj) , one obtains 



EW (F) a £?Ctf) 



m e Aj(F) mhAj (F) 



2M ex j 2 w 2M ex j 2 Wj 



(36) 



where the electron and hole subband shifts are written 
separately. This, in view of Eq. (f3"3"|) . yields the first band 
gap field dependence in the form 



E W(F)kEW 



l-^A^F) 



(37) 



The bang gap decrease with the field in Eq. (|3"T|) is 
stronger than the opposite effect in the negative exciton 
binding energy given (to the same order approximation 
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in field) by Eq. (j35|) . Thus, the first exciton excitation 
energy ([7]) will be gradually decreasing as the perpendicu- 
lar field increases, shifting the exciton absorption peak to 
the red. This is the basic feature of the quantum confined 
Stark effect observed previously in semiconductor nano- 
materials [lO, [Tl|, [73, [73| ■ The field dependences of the 
higher interband transitions exciton excitation energies 
are suppressed by the rapidly (quadratically) increasing 
azimuthal quantization numbers in the denominators of 
Eqs. (P | and 

Lastly, the perpendicular field dependence of the in- 
terband plasmon resonances can be obtained from the 
frequency dependence of the axial surface conductivity 
due to excitons (see Ref. [53[ and refs. therein). One has 



-ifkofj 



(38) 



where and r are the exciton oscillator strength and re- 
laxation time, respectively. The plasmon frequencies are 
those at which the function Re[l/a|f(w)] has maxima. 



Testing it for maximum in the domain < hto < E, 
one finds the first interband plasmon resonance energy to 
be (in the limit r— >oo) 



,(22) 



= 

p 



[E 



(11)12 

exc 



[E 



(22)i 2 

exc 



(39) 



Using the field dependent El xc ' given by Eqs. 0, (j35)) 
and (|3T|) . and neglecting the field dependence of E^xc , 
one obtains to the second order approximation in the field 



E^\F) 



£(11) 
P 



1 +E { g 11} /2E, 
1 +Ei x 2 c/E, 



Kii) 



,.(ii) 

" X C 



Ai(F) 



(40) 



Figure [5] shows the results of our calculations of the 
field dependences for the first bright exciton parame- 
ters in the (11,0) and (10,0) CNs. The energy is mea- 
sured from the top of the first unperturbed hole sub- 
band (as shown in Fig. [2j right panel) . The binding en- 
ergy field dependence was calculated numerically from 
Eq. (|3"4")1 as described above [shown in Fig. [5] (a)]. The 
band gap field dependence and the plasmon energy field 
dependence were calculated from Eqs. (f3"rJ)) and ([4*0]) . re- 
spectively. The zero-field excitation energies and zero- 
field binding energies were taken to be those reported in 
Ref. [HI and in Ref. [Hj], respectively, and we used the 
diameter- and chirality-dependent electron and hole ef- 
fective masses from Ref. [77] . As is seen in Fig.[B](a) and 
(b) , the exciton excitation energy and the interband plas- 
mon energy experience red shift in both nanotubes as 
the field increases. However, the excitation energy red 
shift is very small (barely seen in the figures) due to 
the negative field dependent contribution from the exci- 
ton binding energy. So, E^{F) and E^ l) {F) approach 
each other as the field increases, thereby bringing the 
total exciton energy ^ in resonance with the surface 
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FIG. 6: (Color online) (a),(b) Calculated dependences of the 
first bright exciton parameters in the (11,0) and (10,0) CNs, 
respectively, on the electrostatic field applied perpendicular 
to the nanotube axis. The dimensionless energy is denned as 
[Energy]/ 2^/0, according to Eq. (|17p . The energy is measured 
from the top of the first unperturbed hole subband. 



plasmon mode due to the non-zero longitudinal kinetic 
energy term at finite temperature [78|. Thus, the elec- 
trostatic field applied perpendicular to the CN axis (the 
quantum confined Stark effect) may be used to tune the 
exciton energy to the nearest interband plasmon reso- 
nance, to put the exciton-surface plasmon interaction 
in small-diameter semiconducting CNs to the strong- 
coupling regime. 



C. The optical absorption 

Here we analyze the longitudinal exciton absorption 
line shape as its energy is tuned to the nearest interband 
surface plasmon resonance. Only longitudinal excitons 
(excited by light polarized along the CN axis) couple to 
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the surface plasmon modes as discussed at the very be- 
ginning of this section (see Ref. [56[ for the perpendicular 
light exciton absorption in CNs). We follow the optical 
absorption/emission lineshape theory developed recently 
for atomically doped CNs [8j. (Obviously, the absorption 
line shape coincides with the emission line shape if the 
monochromatic incident light beam is used in the absorp- 
tion experiment.) When the /-internal state exciton is 
excited and the nanotube's surface EM field subsystem 
is in vacuum state, the time-dependent wave function of 
the whole system " exciton+field" is of the form 74 1 



\m) = E c, /( k >*) e " ifi/(k)t/ft i{ 1 /( k )})-i(°}) ( 41 ) 



£ / d W C(k, W ,0e- Mt |{0}) e ,|{l(k, W )}). 



Here, |{l/(k)}) eK is the excited single-quantum Fock 
state with one exciton and |{l(k, to)}} is that with one 
surface photon. The vacuum states are |{0}) ea ; and |{0}) 
for the exciton subsystem and field subsystem, respec- 
tively. The coefficients C/(k, t) and C(k,ui,t) stand for 
the population probability amplitudes of the respective 
states of the whole system. The exciton energy is of the 
form Ef(k) = Ef(k)-ih/T with £>(k) given by Eq. © 
and r being the phcnomcnological exciton relaxation 
time constant [assumed to be such that li/r<£/(k)] to 
account for other possible exciton relaxation processes. 
From the literature we have r p h ~ 30 — 100 fs for the 
exciton-phonon scatterin g l3ll. T<j ~ 50 ps for the exci- 
ton scattering by defects [23, [261, and T ra d ~ 10 ps— 10 ns 
for the radiative decay of excitons [HI]. Thus, the scat- 
tering by phonons is the most likely exciton relaxation 
mechanism. 

We transform the total Hamiltonian (fTj)- (fT0f to the k- 
representation using Eqs. ([5]) and © (see Appendix |A"|) . 
and apply it to the wave function in Eq. (|4"Tj) . We ob- 
tain the following set of the two simultaneous differential 
equations for the coefficients C/(k, t) and C(k, ui, t) from 
the time dependent Schrodinger equation 



complete decay of all initially excited excitons, 

= |C(k,c,^co)| 2 = l^|g( +) (k,k, W )| 2 

/ 

dt' C f (k,t')e-^ E f ( V- h ^ t '/ h . (44) 



Here, the second equation is obtained by the formal inte- 
gration of Eq. (|43|) over time under the initial condition 
C (k, u), 0) = 0. The emission intensity distribution is thus 
related to the exciton population probability amplitude 
C/(k, t) to be found from Eq. (PJ. 

The set of simultaneous equations ([42]) and ([43]) [and 
Eq. (I44|) . respectively] contains no approximations except 
the (commonly used) neglect of many-particle excitations 
in the wave function (|4ip . We now apply these equations 
to the exciton-surface-plasmon system in small-diameter 
semiconducting CNs. The interaction matrix element in 
Eqs. (j4"2"|) and (|4"3")) is then given by the k-transform of 
Eq. (fT3|) , and has the following property (Appendix [C]) 



1 



270 n 



?^ +) (k,k,c)| 2 



1 _. 

— T J (x)p{x) 
Zir 



(45) 



with Tq(x) and p(x) given by Eqs. (flE|) and (TIT))) , respec- 
tively. We further substitute the result of the formal inte- 
gration of Eq. (43]) [with C(k,w,0) = 0] into Eq. [|l2]). use 
Eq. (j45]l with p(x) approximated by the Lorentzian (|21jl , 
calculate the integral over frequency analytically, and dif- 
ferentiate the result over time to obtain the following sec- 
ond order ordinary differential equation for the exciton 
probability amplitude [dimcnsionless variables, Eq. (|17|) ] 

C f ([3) + lAx p -Ae f +i(x p -e f )}C f (f3) + (Xf/2) 2 Cf(f3) = 0, 

where X f = [2Ax p T f(x p )] 1 ^ 2 with f f(x p ) =f l(x p )p(x p ), 
Aef = h/2-f T, P = 2^ t/h is the dimensionless time, and 
the k-dependence is omitted for brevity. When the total 
exciton energy is close to a plasmon resonance, Ef~x p , 
the solution of this equation is easily found to be 



dugf\k,k'^) (42) 



x C(k',uj,t)e- iwt 



<7(k>,t)e-- t 5 kk , = - i£ [g( +) (k,k»]* (43) 

/ 

x CXk,t)e-^W*/ n . 



The ^-symbol on the left in Eq. (j43|) ensures that the 
momentum conservation is fulfilled in the exciton-photon 
transitions, so that the annihilating exciton creates the 
surface photon with the same momentum and vice versa. 
In terms of the probability amplitudes above, the exci- 
ton emission intensity distribution is given by the final 
state probability at very long times corresponding to the 



2 1 /Sx 2 - X 2 t 



f , 



8x 



6x2 ~ X f, 



-(Sx+y/Sx 2 -X^l3/2 



(46) 



where Sx = Ax n — As f > and X f = [2Ax p T f (e f )]^ 2 . 
This solution is valid when Sf ~ x p regardless of the 
strength of the exciton-surface-plasmon coupling. It 
yields the exponential decay of the excitons into plas- 
mons, \Cf((3)\ 2 w exp[— f /(e/)/3], in the weak coupling 
regime where the coupling parameter (Xf/Sx) 2 <C 1. 
If, on the other hand, (Xf/5x) 2 ^ 1, then the strong 
coupling regime occurs, and the decay of the excitons 
into plasmons proceeds via damped Rabi oscillations, 
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|C/(/3)| 2 «exp(-fe/3)cos 2 (X / /3/2). This is very similar 
to what was earlier reported for an excited two-level atom 
near the nanotube surface [lO, IH, [45[ . Note, however, 
that here we have the exciton-phonon scattering as well, 
which facilitates the strong exciton-plasmon coupling by 
decreasing Sx in the coupling parameter. In other words, 
the phonon scattering broadens the (longitudinal) exci- 
ton momentum distribution [81] , thus effectively increas- 
ing the fraction of the excitons with Ef~x p . 

In view of Eqs. (|43)l and (146[) . the exciton emission 
intensity (|44[) in the vicinity of the plasmon resonance 
takes the following (dimensionless) form 



I(x) « T (e f ) 



f 



dpC f (f3)e 



i(x—6f-\-iAef)0 



(47) 



where I(x) = 2jqI(ui)/?i and Iq = T t(Ef)/2%. After some 
algebra, this results in 



I(x) 



I {e f )[(x-e f f 



[(x - e f ) 2 - X 2 /4} 2 + (x - e f ) 2 (Ax 2 p + As 2 ) ' 



(48) 

where Ax 2 > Ae 2 . The summation sign over the exciton 
internal states is omitted since only one internal state 
contributes to the emission intensity in the vicinity of 
the sharp plasmon resonance. 

The line shape in Eq. (|48p is mainly determined by the 
coupling parameter (Xf/Ax p ) 2 . It is clearly seen to be 
of a symmetric two-peak structure in the strong coupling 
regime where (X f / Ax p ) 2 ^> 1. Testing it for extremum, 
we obtain the peak frequencies to be 



Xi, 2 = Ef ± 




Axj_ 
X f 



[terms ~ (Ax p ) 2 (Ae f) 2 / Xj are neglected], with the Rabi 
splitting x\ — X2 ~ Xf . In the weak coupling regime where 
(Xf/Ax p ) 2 <C 1, the frequencies x\ and xi become com- 
plex, indicating that there are no longer peaks at these 
frequencies. As this takes place, Eq. (|4"5|) is approximated 
with the weak coupling condition, the fact that x~£f, 
and X 2 = 2Ax p T f(e f), to yield the Lorentzian 



I (e f )/[l + (Ae f /Ax p ) 2 



(x-e f ) 2 + Tf(ef)/2J1 + (Ae f /Ax p ) 



I(x) 



peaked at x = Ef, whose half- width- at-half-maximum is 
slightly narrower, however, than Tf(ef)/2 it should be if 
the exciton-plasmon relaxation were the only relaxation 
mechanism in the system. The reason is the competing 
phonon scattering takes excitons out of resonance with 
plasmons, thus decreasing the exciton-plasmon relax- 
ation rate. We therefore conclude that the phonon scat- 
tering does not affect the exciton emission/ absorption 
line shape when the exciton-plasmon coupling is strong 



(it facilitates the strong coupling regime to occur, 
however, as was noticed above), and it narrows the 
(Lorentzian) emission/absorption line when the exciton- 
plasmon coupling is weak. 

Calculated exciton emission/absorption lincshapes, as 
given by Eq. (|48|) for the CNs under consideration, are 
shown in Fig. [7] (a) and (b). The exciton energies are 
assumed to be tuned, e. g. by means of the quantum 
confined Stark effect discussed in Sec. IIIII B. to the near- 
est plasmon resonances (shown by the vertical dashed 
lines in the figure). We used r p h = 30 fs as reported 
in Ref. [23]. The line (Rabi) splitting effect is seen to be 
~ 0.1 eV, indicating the strong exciton-plasmon coupling 
with the formation of the mixed surface plasmon-exciton 
excitations. The splitting is larger in the smaller diame- 
ter nanotubes, and is not masked by the exciton-phonon 
scattering. 



IV. CONCLUSIONS 

We have shown that the strong exciton-surface- 
plasmon coupling effect with characteristic exciton ab- 
sorption line (Rabi) splitting ~ 0.1 eV exists in small- 
diameter (< 1 nm) semiconducting CNs. The splitting 
is almost as large as the typical exciton bindin g e n- 
ergies in such CNs (~ 0.3 - 0.8 eV [H E EI Hi), 
and of the same order of magnitude as the exciton- 
plasmon Rabi splitting in organic semiconductors (~ 
180 meV 37]). It is much larger than the exciton- 
polariton Rabi splitting in semiconductor microcavities 
(~ 140 - 400/zeV [H il HI), or the exciton-plasmon 
Rabi splitting in hybrid semiconductor-metal nanoparti- 
cle molecules [38| . 

Since the formation of the strongly coupled mixed 
exciton-plasmon excitations is only possible if the ex- 
citon total energy is in resonance with the energy of 
an interband surface plasmon mode, we have analyzed 
possible ways to tune the exciton energy to the near- 
est surface plasmon resonance. Specifically, the exciton 
energy may be tuned to the nearest plasmon resonance 
in ways used for the excitons in semiconductor quantum 
microcavities — thermally (by elevating sample tempera- 
ture) [67], [H, [69[ , and/or electrostatically [70. l7ll. WX [73j 
(via the quantum confined Stark effect with an exter- 
nal electrostatic field applied perpendicular to the CN 
axis) . The two possibilities influence the different degrees 
of freedom of the quasi-lD exciton — the (longitudinal) 
kinetic energy and the excitation energy, respectively. 

We have studied how the perpendicular electrostatic 
field affects the exciton excitation energy and interband 
plasmon resonance energy (the quantum confined Stark 
effect). Both of them are shown to shift to the red due 
to the decrease in the CN band gap as the field increases. 
However, the exciton red shift is much less than the plas- 
mon one because of the decrease in the absolute value 
of the negative binding energy, which contributes largely 
to the exciton excitation energy. The exciton excitation 
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FIG. 7: (Color online) (a),(b) Exciton absorption/emission 
lineshapes as the exciton energies are tuned to the nearest 
plasmon resonance energies (vertical dashed lines in here; see 
Fig. |3]and left panels in Fig.g} in the (11,0) and (10,0) nano- 
tubes, respectively. The dimensionless energy is defined as 
[Energy]/ 2^0, according to Eq. (fTT)l . 



energy and interband plasmon energy approach as the 
field increases, thereby bringing the total exciton energy 
in resonance with the plasmon mode due to the non-zero 
longitudinal kinetic energy term at finite temperature. 

Lastly, the noteworthy message we would like to de- 
liver in this paper is that the strong exciton-surface- 
plasmon coupling we predict here occurs in an individ- 
ual CN as opposed to various artificially fabricated hy- 
brid plasmonic nanostructures mentioned above. We 
strongly believe this phenomenon, along with its tun- 
ability feature via the quantum confined Stark effect we 
have demonstrated, opens up new paths for the develop- 
ment of CN based tunable optoelectronic device applica- 
tions in areas such as nanophotonics, nanoplasmonics, 
and cavity QED. One straightforward application like 



this is the CN photoluminescence control by means of 
the exciton-plasmon coupling tuned electrostatically via 
the quantum confined Stark effect. This complements 
the microcavity controlled CN infrared emitter appli- 
cation reported recently [25j], offering the advantage of 
less stringent fabrication requirements at the same time 
since the planar photonic microcavity is no longer re- 
quired. Electrostatically controlled coupling of two spa- 
tially separated (weakly localized) excitons to the same 
nanotube's plasmon resonance would result in their en- 
tanglement [HI [H| , the phenomenon that paves the 
way for CN based solid-state quantum information ap- 
plications. Moreover, CNs combine advantages such as 
electrical conductivity, chemical stability, and high sur- 
face area that make them excellent potential candidates 
for a variety of more practical applications, including ef- 
ficient solar energy conversio n [7| , energy storage [12j, 
and optical nanobiosensorics [861 ] - However, the photo- 
luminescence quantum yield of individual CNs is rela- 
tively low, and this hinders their uses in the aforemen- 
tioned applications. CN bundles and films are proposed 
to be used to surpass the poor performance of individ- 
ual tubes. The theory of the exciton-plasmon coupling 
we have developed here, being extended to include the 
inter-tube interaction, complements currently available 
'weak-coupling' theories of the exciton-plasmon interac- 
tions in low-dimensional nanostructures [38, 87] with the 
very important case of the strong coupling regime. Such 
an extended theory (subject of our future publication) 
will lay the foundation for understanding inter-tube en- 
ergy transfer mechanisms that affect the efficiency of op- 
toelectronic devices made of CN bundles and films, as 
well as it will shed more light on the recent photolu- 
minescence experiments with CN bundles (88l . l89l | and 
multi- walled CNs [9(J, revealing their potentialities for 
the development of high-yield, high-performance opto- 
electronics applications with CNs. 
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APPENDIX A: EXCITON INTERACTION WITH 
THE SURFACE ELECTROMAGNETIC FIELD 

We follow our recently developed QED formalism to 
describe vacuum-type EM effects in the presence of quasi- 
1D absorbing and dispersive bodies El, E|, IH, HI, 
EI- The treatment begins with the most general EM 
interaction of the surface charge fluctuations with the 
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quantized surface EM field of a single- walled CN. No ex- 
ternal field is assumed to be applied. The CN is modelled 
by a neutral, infinitely long, infinitely thin, anisotropi- 
cally conducting cylinder. Only the axial conductivity of 
the CN, a zz , is taken into account, whereas the azimuthal 
one, er w , is neglected being stron g ly suppressed b y the 



transverse depolarization effect [49|, [50|, iMl 152|, [53J, [5j 
Since the problem has the cylindrical symmetry, the or- 
thonormal cylindrical basis {e r ,e v ,e z } is used with the 
vector e z directed along the nanotube axis as shown in 
Fig. [1] The interaction has the following form (Gaussian 
system of units) 



Hi, 



H. 



n,z 



(1) 



H 



(2) 



(Al) 



pW_|A(n + fW) 



where c is the speed of light, m,, g;, f„ , and are, re- 
spectively, the masses, charges, coordinate operators and 
momenta operators of the particles (electrons and nu- 
cleus) residing at the lattice site n = R„ = {Rcn , 'fin, z n} 
associated with a carbon atom (see Fig. [TJ on the sur- 
face of the CN of radius Rc n ■ The summation is taken 
over the lattice sites, and may be substituted with the 
integration over the CN surface using Eq. ([3]). The vec- 
tor potential operator A and the scalar potential opera- 
tor (p represent the nanotube's transversely polarized and 
longitudinally polarized surface EM modes, respectively. 
They are written in the Schrodinger picture as follows 



A(n) = / du — E (n, u) + h.c, 



(A2) 



-V„<£(n)= / rfcjE (n,Lu) + h.c. (A3) 
Jo 

We use the Coulomb gauge whereby V n -A(n) = 0, or 

equivalently, \pn , A(n + r$)] = 0. 

The total electric field operator of the CN-modified EM 
field is given for an arbitrary r in the Schrodinger picture 
by 



E(r) = / dwE(r,w) + h.c. 



(A4) 



doj\E (r, co) + E (r, w)] + h.c. 



with the transversely (longitudinally) polarized Fourier- 
domain field components defined as 

E ±(ll) (r,cj) =fdr'6 ± W(r-r')-E(r',u>), (A5) 



(A6) 



5^{v)=8 a p8{r)-5^{v) 



are the longitudinal and transverse dyadic ^-functions, 
respectively. The total field operator (|A4p satisfies the 
set of the Fourier-domain Maxwell equations 



V x E(r,w) = ikR(T,u), 



(A7) 



47T 



V x H(r, lu) = -ik E(r, u) + — I(r, w), (A8) 

c 

where H = (ifc) _1 V x E is the magnetic field operator, 
k = uj/c, and 



(A9) 



is the exterior current operator with the current density 
defined as follows 



where 



* hujRea zz (R C N,u) t, x , Ain >> 
J(n,w) = y f(n,uj)e z (A10) 

to ensure preservation of the fundamental QED equal- 
time commutation relations (see, e.g., (46|) for the EM 
field components in the presence of a CN. Here, a zz is 
the CN surface axial conductivity per unit length, and 
/(n, w) along with its counter-part p (n, u>) are the scalar 
bosonic field operators which annihilate and create, re- 
spectively, single-quantum EM field excitations of fre- 
quency u> at the lattice site n of the CN surface. They 
satisfy the standard bosonic commutation relations 

[/(n 1 w),/t(m ) wO]=inm«(w-w'). (All) 
[/(n,w),/(m,</)] = [/t( n , w ),/t(m,^)] =0. 
One further obtains from Eqs. (jA"7|) - (|A10p that 

4-7T 

E(r,w) =ik— ^G(r,n,w)-J(n,w), (A12) 

n 

and, according to Eqs. (|A4j) and (|A5|l . 

E ±(ll) (r,w) =ik— V ^G(r,n,u)-3(n,Lj), (A13) 
c ^— ' 

n 

where "'"G and "G are the transverse part and the lon- 
gitudinal part, respectively, of the total Green tensor 
G = ^G + "G of the classical EM field in the presence 
of the CN. This tensor satisfies the equation 

J2 (V x V x - k 2 ) za G az (r, n, w) = S(r - n) (A14) 

a—r,ip,z 

together with the radiation conditions at infinity and the 
boundary conditions on the CN surface. 
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All the 'discrete' quantities in Eqs. (|A"9 |) - (|A14[) may 
be equivalently rewritten in continuous variables in view 
of Eq. ([3]). Being applied to the identity 1 = J2 m ^nm, 
Eq. © yields 



This requires to redefine 



(A15) 



/(n,w) = v^SbARn.w), /t(n,w) = ^P(R n ,uj) 

(A16) 

in the commutation relations (|A11|) . Similarly from 
Eq. (|AT2]1 . in view of Eqs. ©, ([ATI))) and (jA"l%|) . one 

obtains 



G(r,n,w) = v%G(r,R n ,w), 



(A17) 



which is also valid for the transverse and longitudinal 
Green tensors in Eq. (|A13[) . 

Next, we make the series expansions of the interactions 
Hint and H^l in Eq. (|A1|) about the lattice site n to the 
first non- vanishing terms, 

9i a f„\ 



rV(l) 
n int 



' — ' Ul-iC 



E 



9i 



2m,c 2 



A 2 (n), (A18) 



and introduce the single-lattice-site Hamiltonian 

H n =e a \0)(0\ + $>o+^/)|/) (f\ 
f 

with the completeness relation 

|o)<o| + £|/)</| = j. 



(A19) 



(A20) 



(A21) 



Here, eo is the energy of the ground state |0) (no exciton 
excited) of the carbon atom associated with the lattice 
site n, eo+htuf is the energy of the excited carbon atom in 
the quantum state |/) with one /-internal-state exciton 
formed of the energy Eill = huf. In view of Eqs. (|A20[) 
and (|A21[) . one has 



dr 



(i) 



- \r M H 1 - — / If (0 if 1 f 

j. — t- L x n j-^nj — . fc [ x n j^nJJ 

at «ra in 



^X>/(<0|*«|/)iW- (f\*$\0)Bl f ) (A22) 



m 

2 



B n j = |0)(/| and B^ , = |/)(0| create and annihilate, 
respectively, the /-internal-state exciton at the lattice 
site n, and exciton-to-exciton transitions are neglected. 
In addition, we also have 



SaS, 



ij 0ol(3 



(A24) 



where a, f3 = r,tp,z. Substituting these into Eqs. (|A18|) 
and (I A19|) [commutator (I A24[) goes into the second term 
of Eq. (j A18|) which is to be pre-transformed as follows 
J2i,j,a,0 Qi<ljM n )aM n )i3SijS a ^/2m i c 2 }, one arrives at the 
following (electric dipole) approximation of Eq. (|A1[) 



Hint = H 



(i) 



r(2) 



H 

int int 



(A25) 



■£^-A(n) 



n.f 



B 



n,f~ 



B»,/+^di-A(n) 



+ ]Td£.V n ^n) (Bl f + B n . 



n,/' 

with d£ = (0|d n |/) = (/|d n |0), where d n = Ei^n « 
the total electric dipole moment operator of the particles 
residing at the lattice site n. 

The Hamiltonian (|A25|1 is seen to describe the vacuum- 
type exciton interaction with the surface EM field (cre- 
ated by the charge fluctuations on the nanotube surface) . 
The last term in the square brackets does not depend on 
the exciton operators, and therefore results in the con- 
stant energy shift which can be safely neglected. We then 
arrive, after using Eqs. (|A2|) . (|A3]l . (|A10|) . and (|A13[) , at 
the following second quantized interaction Hamiltonian 



/>oo 

H mt = V / ^[g^ +) (n,m, W )Bt 
Jo 



g ( r } (n, m, Lo)B aJ } /(m, w) + h.c, (A26) 



where 



gf ] (n, m, w) = g j (n, m, w) ± — g'i (n, m, cj) (A27) 



with 



g| (ll) (n, m,w) = Trftw Re a zz (R C N, w) 

x E(d0a i(ll) G a ,(n,m, w ), (A28) 



a— r,(/?,2 



and 



and 



53(<0|?«|/)B n , / + (/|rW|0)Bi i/ V ±(ll) G QZ (n,m,w) =JdrS^\n-r) G^(r,m, £ 



(A23) 

where (0|fn' ) |/) = (/|r^'|0) in view of the hermitian and 
real character of the coordinate operator. The operators 



(A29) 

This yields Eqs. (fl"U|) - (fT!?]) after the strong transverse de- 
polarization effect in CNs is taken into account whereby 



if 



(<0 
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APPENDIX B: GREEN TENSOR OF THE 
SURFACE ELECTROMAGNETIC FIELD 

Within the model of an infinitely thin, infinitely long, 
anisotropically conducting cylinder we utilize here, the 
classical EM field Green tensor is found by expanding the 
solution to the Green equation (|A14|) in series in cylin- 
drical coordinates, and then imposing the appropriately 
chosen boundary conditions on the CN surface to de- 
termine the Wronskian normalization constant (see, e.g., 
Ref. |3). 

After the EM field is divided into the transversely 
and longitudinally polarized components according to 
Eqs. (|X3] l -(|A6] l . the Green equation (|A14|) takes the form 



(Vx Vx -fc 2 ) 2Q [ ± G QZ (r,n,w)+llG Q2 (r,n,w) 

a— r,tp,z 

= S(r-n) (Bl) 
with the two additional constraints, 

X) V«, J -G a ,(r,n,w) = (B2) 



and 



e a /3-yV/3 "G 7Z (r, n,w) =0, 



(B3) 



where e Q/ 3 7 is the totally antisymmetric unit tensor of 
rank 3. Equations (|B2|) and (|B3[) originate from the 
divergence-less character (Coulomb gauge) of the trans- 
verse EM component and the curl-less character of the 
longitudinal EM component, respectively. The trans- 
verse ± G az and longitudinal "G Q2 Green tensor compo- 
nents are defined by Eq. (|A29j) which is the corollary of 
Eq. using the Eqs. (|AT2|) and (|Al3]l . Equation (|BT]) 
is further rewritten in view of Eqs. (|B2|) and (|B3[) . to give 
the following two independent equations for G zz and 
"G Z2 we need 

(A + fc 2 ) ^G zz (r, n, w) = -<£(r - n) , (B4) 

k 2 WG zz {r,n,uj) =-4(r-n) (B5) 

with the transverse and longitudinal delta-functions de- 
fined by Eq. ([Ae]) , 

We use the differential representations for the trans- 
verse ± G ZZ and longitudinal "G 22 Green functions of the 
following form [consistent with Eq. (|A29[) ] 

± G zz (r,n,u;) = f-Lv»V, + lW,n,w), (B6) 



HG„(r,n,w) = -^jV,V»s(r,n,w), (B7) 

where <?(r, n, cj) is the scalar Green function of the 
Helmholtz equation (|B4[) . satisfying the radiation con- 
dition at infinity and the finiteness condition on the axis 



of the cylinder. Such a function is known to be given by 
the following series expansion 



ff(r,n,w) 



/So e ife l r ~ R "l 



So 



4tt 



(2tt 



\2 



„ip{lp—iPn) 



x / dhI p {vr)K p {vR CN )e lh ^ z ~ z -\ r < R CN , 



(B8) 



where I p and K p arc the modified cylindric Bessel func- 
tions, v = v(h, uj) = \Jh 2 ~ fc 2 , and we used the property 
(IA17P to go from the discrete variable n to the corre- 
sponding continuous variable. The integration contour 
G goes along the real axis of the complex plane and en- 
velopes the branch points ±fc of the integrand from below 
and from above, respectively. For r > Rcn> the function 
g(r, n, u>) is obtained from Eq. (|B8[) by means of a simple 
symbol replacement I p <-* K p in the integrand. 

The scalar function (|B8[) is to be imposed the bound- 
ary conditions on the CN surface. To derive them, we 
represent the classical electric and magnetic field compo- 
nents in terms of the EM field Green tensor as follows 



E a {Y,uS) = ik- L G az {r,n,w), 



(B9) 



H a (r,w) = ~ J2 £ ^V^ 7 (^)- (BIO) 

These are valid for r =Z n under the Coulomb-gauge con- 
dition. The boundary conditions arc then obtained from 
the standard requirements that the tangential electric 
field components be continuous across the surface, and 
the tangential magnetic field components be discontinu- 
ous by an amount proportional to the free surface current 
density, which we approximate here by the (strongest) 
axial component, <j zz {Rcn of the nanotube's surface 
conductivity. Under this approximation, one has 



E z \+-E z \_=E v \ + -E v \_ = 0, 

H z \+-H z \- = 0, 
Air 

HJ+-MJ- = —<r zz {u;)E z \ RGN , 



(Bll) 
(B12) 

(B13) 



where ± stand for r = Rqn i £ with the positive in- 
finitesimal e. In view of Eqs. (|B9|) . (|B10j) and (|B6|) . the 
boundary conditions above result in the following two 
boundary conditions for the function (|B8j) 



9\ 



Or 



dg 
Or 



0, 

4nia zz (uj) ( 9_ 
dz 



(B14) 



k z \g\ RcN . (B15) 



We see that Eq. (|Bl"4)) is satisfied identically. Eq. (|B"15|) 
yields the Wronskian of modified Bessel functions on the 
left, W[I P (x) , K p (x)} = I p {x)K' p (x) - K p (x)I' p {x) = -l/x, 
which brings us to the equation 

_ 1 = AniaM v 2j p{vRcN)Kp{vRcN y (m6) 
Rcn U) 



16 



This is nothing but the dispersion relation which deter- 
mines the radial wave numbers, h, of the CN surface 
EM modes with given p and ui. Since we are interested 
here in the EM field Green tensor on the CN surface 
[see Eq. (|A28|) ] . not in particular surface EM modes, 
we substitute I p (vRcn)K p (vRcn) from Eq. (|B16jl into 
Eq. (|B8p with r ~ Rcn- This allows us to obtain the 
scalar Green function of interest with the boundary con- 
ditions (|B14p and (|B15|) taken into account. We have 



g(R,n,w) = - 



iuJy/So6(ip - (fin) 

8ir 2 a zz (uj)RcN 



dh- 



^ih(z — z n ) 

k 2 -h 2 



(B17) 



where R = {Rcn, f, z} is an arbitrary point of the cylin- 
drical surface. Using further the residue theorem to cal- 
culate the contour integral, we arrive at the final expres- 
sion of the form 



(B18) 



8ira zz {LJ)RcN 



which yields 



x G zz {R,n,u) = 0, 
"G ZZ (R, n,w) = g(R,n,w), 



(B19) 
(B20) 



in view of Eqs. ()B6} and (fB7j. 

The fact that the transverse Green function (|B19p 
identically equals zero on the CN surface is related to the 
absence of the skin layer in the model of the infinitely thin 
cylinder (see, e.g., Ref. [82]). In this model, the trans- 
verse Green function is only non-zero in the near-surface 
area where the exciton wave function goes to zero. Thus, 
only longitudinally polarized EM modes with the Green 
function (|B20|) contribute to the exciton surface EM field 
interaction on the nanotube surface. 



APPENDIX C: DIAGONALIZATION OF THE 
HAMILTONIAN ©-(T13]) 

We start with the transformation of the total Hamil- 
tonian (fT|)- p3|) to the k-representation using Eqs. (0) 
and ((9]). The unperturbed part presents no difficulties. 
Special care should be given to the interaction matrix el- 
ement g^(n, m, uj) in Eq. (fl~3|) . In view of Eqs. (|B20|) . 
(|B18P and ([3]), one has explicitly 

gf^k'.u,) = 1 ^g^(n,m,w) e~ lkn+lk ' m (CI) 



_ ± Wy/7r7kjReo- zz (^) ^{y^T 
2ttc(t zz (uj)Rcn N 



x / dz n dz m e 



iuj\z n — z m \j c—ik z z n +ik' z z n 



where we have also taken into account the fact that the 
dipole matrix element (d£) z = (0|(d n ) z |/) is actually the 
same for all the lattice sites on the CN surface in view of 
their equivalence. As a consequence, (d£) z = d{/N with 

d£=En<0|(d n ),|/>. . . 

The integral over ip in Eq. (|Cip is taken in a standard 
way to yield 

d Vn d^ m 6(tp n - ^ m )e- ik ^ +tk >™ =2tt4 w . (C2) 



The integration over z is performed by first writing the 
integral in the form 



L — *oo 



L/2 f L/2 



-L/2 J-L/2 



(L being the CN length), then dividing it into two parts 
by means of the equation 

e H*»-*m|/c = 6»(z„-z m )e^ (z "~ Zm)/c 
+ 6{z m -z n )e- l ^ z -- z ^l c 1 

and finally by taking simple exponential integrals with 
allowance made for the formula 



lim 



2sin[L(fc z - k' z )/2] 



l^oo L(k z — k' z ) 
After some simple algebra we obtain the result 

'\z n — z m \/c— ik z z n +ik' z z m 



dz n dz m 6 



(C3) 



2iuj /c 



lim L 2 1 1 — 

l^oo I L[k 2 - (w/c) 2 l 



5k. 



In view of Eqs. fC2[) and ((C3|) . the function (jClj) takes 
the form 



gf J (k,k» = ±- 



x lim < 1 — 



' d{ y irSohuj Re <r zz (i 
{2n) 2 ca zz (uj)RcN 

2iio I c 



(C4) 



i^oo (_ L [k 2 — (w/c) 2 



<5kk'- 



We have taken into account here that Sk^k' 5k z k' z = <5kk', 
as well as the fact that (R C nL/NS ) 2 = 1/(2tt) 2 . This 
can be further simplified by noticing that only absolute 
value squared of the interaction matrix element matters 
in calculations of observables. We then have 



1 



2iuj/c 



L[k 2 -(co/c) 2 



1 



a 

' 2 



1 



n 



u l + or 



with u = (ck z /uj) 2 — 1, and a = (2c/Loj) 2 being the small 
parameter which tends to zero as L — * oo. Using further 
the formula (see, e.g., Ref. [6Cl |) 



5(u) 



1 



lim 



a 



7T a—>o vr + or 
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and the basic properties of the (5-function, we arrive at 

2 17 1 



lim 

L — >oc 



1 



2iuj/c 



L[k2-(u/cY 



= 1 



(C5) 



x [S(uj + ck z ) + S(w — ck x )] 



We also have 



y/Recr zz (LLi) 



= Re 



o zz {uS) 



(C6) 



Equation (|C4[) . in view of Eqs. (|C5[) and (|C6[) . is rewritten 
effectively as follows 



with 



x\ 1 + 



' d{ y/ TrSphuj Re[l/a zz (uj)} 
(2tt) 2 cR C n 

-!■ [8{uj + ck z ) + 5{uj-ck z )] 



(C7) 



(C8) 



In terms of the simplified interaction matrix element 
(|C7[) , the k-representation of the Hamiltonian (fljl- fT^)) 
takes the following (symmetrized) form 



(C9) 



where 



H k = E f (k) [Bl f B kJ + B\ f B^ f ) (CIO) 

/ 

-Jluhtu [/t(k,w)/(k,w) + /t(-k,w)/(-k, 



h.c. 



with Df(uj) given by Eq. (|C8[) . To diagonalize this Hamil- 
tonian, we follow Bogoliubov's canonical transformation 
technique (see, e.g., Ref. |60(). The canonical transfor- 
mation of the exciton and photon operators is of the form 



B kJ = [^(k,o;/)i(k) + » p (k,w / )#(-k)] , (Cll) 
/(k,u/) = J] [^(k, W )4(k) +<(k, w )|t(_ k )] , (C12) 

where the new operators, £ p (k) and £^(k) = [^(k)]^, 
annihilate and create, respectively, the coupled exciton- 
photon excitations of branch fi on the nanotube surface. 
They satisfy the bosonic commutation relations of the 
form 



4(k),£t(k') 



= <?W"5kk', 



(C13) 



which, along with the reversibility requirement of 
Eqs. (jClip and (|C12[) . impose the following constraints 
on the transformation functions and v M 

53 KO^/W'O 4 '"/) -^( k > w /K'( k > w /)] 
/ 



+ / duj [u p (k, u)u*, (k, oj) - v* (k, u)v^ (k, u>)] = 6^ , 
Jo 

51 K( k ' w /) u M( k . w /') -^( k ^/)^( k ' w /')] = 

53 [ u ;(k, u ) Uft (k, lu') - ^(k, W )«„(k, «/)] = «(w - w')- 



Here, the first equation guarantees the fulfilment of the 
commutation relations (|C13|) . whereas the second and 
the third ensure that Eqs. (jCllj) and (|C12|) are inverted 
to yield £ M (k) as gi ven by Eq. (fT5)) . Other possible com- 
binations of the transformation functions are identically 
equal to zero. 

The proper transformation functions that diagonalize 
the Hamiltonian (|C10[) to bring it to the form ([14]) , are 
determined by the identity 



7iuv(k)4(k)= 4(k),£ k 



(C14) 



Putting Eqs. §T5§ and (IC10|) into Eq. (1C14|) and using 
the bosonic commutation relations for the exciton and 
photon operators on the right, one obtains (k-argument 
is omitted for brevity) 



(hujf, - E f )u*(cjf) =-i duD f ((j) [u^uj) - v*(w)] , 
Jo 

/>oo 

(hu^+ E f )v^{Lu f ) = i / duD f (u) [u M H - v*(w)] , 
Jo 

h (u^ - uj) u^u) = ^53 D f(u) [ u *^f) + v »(uf)] > 
/ 

h (w p + ui) u*(w) = i53 %H Kl( w /) + v ^f)] ■ 
f 

These simultaneous equations define the complex trans- 
formation functions and v p uniquely. They also define 
the dispersion relation (the energies tku^, /j, = 1, 2) of the 
coupled exciton-photon (or exciton-plasmon, to be exact) 
excitations on the nanotube surface. Substituting it M and 
v* from the third and forth equations into the first one, 
one has 



AE 



f 



u\D f (u)\ 



2 1 



Uuj^ + Ef J ?i(lj2 - w 2) 



u*Jw f ) = 0, 



whereby, since the functions u* are non-zero, the disper- 
sion relation we are interested in becomes 



(ft^) 2 



Jo Ho% ~ w 2 ) 



(C15) 
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The energy Eq of the ground state of the coupled exciton- 
plasmon excitations is found by plugging Eq. (fT5)) into 
Eq. (fl4|) and comparing the result with Eqs. (|C9j) and 
(|UT0| . This yields 



E = -^2^(k) 



k, n=l,2 



El^(k^/)| 2 



<2w |i» M (k, cl>)|" 



Using further Df{u>) as explicitly given by Eq. (|C8 
the dispersion relation (|C15|) is rewritten as follows 



£/fr>l4l 2 



(III; 



w 4 Re[l/a zz (a;)] 



7r(c|fc z |) 5 Re[l/a z ,(c|fc z 
ul-(c\k z \) 2 



Here we have taken into account the general property 
Czz(^) = o'tzi - UJ )^ which originates from the time- 
reversal symmetry requirement, in the second term on 
the right hand side. This term comes from the two delta 
functions in \Df(w)\ 2 , and describes the contribution of 
the spatial dispersion (wave-vector dependence) to the 
formation of the exciton-plasmons. We neglect this term 
in what follows because the spatial dispersion is neglected 
in the nanotube's axial surface conductivity in our model, 
and, secondly, because it is seen to be very small for not 
too large excitonic wave vectors. Thus, converting to the 
dimensionless variables (|17p . we arrive at the dispersion 
relation (| 16[) with the exciton spontaneous decay (recom- 
bination) rate and the plasmon DOS given by Eqs. (fl8|) 
and (fl9|) . respectively. 

Lastly, bearing in mind that the delta functions in 
|-D/(w)| 2 are responsible for the spatial dispersion which 
we neglect in our model, and therefore dropping them 
out from the squared interaction matrix element (|C7[) . 
we arrive at the property (|45|) . 



APPENDIX D: EFFECTIVE LONGITUDINAL 

POTENTIAL IN THE PRESENCE OF THE 
PERPENDICULAR ELECTROSTATIC FIELD 



Here we analyze the set of equations (|26|) - (l28l) . and 
show that the attractive cusp-type cutoff potential ([HI]) 
with the field dependent cutoff parameter ([52]) is a uni- 
formly valid approximation for the effective electron- hole 



Coulomb interaction potential ([29)) in the exciton binding 
energy equation (f!?8"|) . 

We rewrite Eqs. (|2l))) and ^Ff\ in the form of a single 
equation as follows 



(Dl) 



— +P COStp ) 1p(tp) = 



and 



Here, ip = tp eyh , ip = ip eyh , q = RcNy/2rn e ,h£e,h/h: 
p = ±2em e ^RcN^ r /^ 2 with the (+)-sign to be taken for 
the electron and the (-)-sign to be taken for the hole. We 
are interested in the solutions to Eq. (jDlj) which satisfy 
the 27r-periodicity condition ip(tp) — ip(tp + 2tt). The 
change of variable tp = 2t transfers this equation to the 
well known Mathieu's equation (see, e.g., Refs. [83l. I&jj ) . 
reducing the solution's period by the factor of two. The 
exact solutions of interest are, therefore, given by the odd 
Mathieu functions se2m+2{t — tp/2) with the cigen values 
&2m+2, where m is a nonnegative integer (notations of 
Ref. [83]). These are the solutions to the Sturm-Liouville 
problem with boundary conditions on functions, not on 
their derivatives. 

It is easier to estimate the z-dependence of the poten- 
tial (|29[) if the functions ip e ,h{}Pe,h) are known explicitly. 
So, we do solve Eq. (|D1|) using the second order pertur- 
bation theory in the external field (the term p cos tp) . The 
second order field corrections are also of practical impor- 
tance in the most of experimental applications. 

The unperturbed problem yields the two linearly inde- 
pendent normalized eigen functions and the eigen values 
as follows 



exp(±ijtp) 



q = j 



Rcn 



(D2) 

with j being a nonnegative integer. The energies s^lij) 

are doubly degenerate with the exception of efl(0) = 0, 
which wc will discard since it results in the zero unper- 
turbed band gap according to Eq. |5]). The perturba- 
tion p cos tp does not lift the degeneracy of the unper- 
turbed states. Therefore, we use the standard nonde- 
generate perturbation theory with the basis wave func- 
tions set above (plus sign selected for defmitcness) to 
calculate the energies and the wave functions to the sec- 
ond order in perturbation. The standard procedure (see, 
e.g., Ref. [H) yields 



Ipj e,h{Ve,h) = ( 1 - 



- 2) 



l e,h 



e 2 R' 



CN 



i(j I) 2 



± 



jT [U 

(<Pe,h) 



2fr 4 



F 



Tpfihi^h) 



(D3) 
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tf(j-2)tf(j 



3 )V , j ( -2e^(^^0 



i e -"-CAT 



[(j-i) 2 -j 2 ][(i-2) 2 -f 



[(j + l) 2 -J 2 ][(j + 2) 2 -f 



Here, j is a positive integer, and the theta-functions en- 
sure that j = 1 is the ground state of the system. The 
corresponding energies are as follows 



Se,h 



h 2 j 2 m eJl e 2 R 4 
2m e j l R 2 CN 2h 



CN 3 p2 



(D4) 



with Wj given by Eq. (|3"3")h thus, according to Eq. ([H]), 
resulting in the nanotube's band gap as given by Eq. (|36|) . 

From Eq. (fD3|) . in view of Eq. (ID2]) . we have the fol- 
lowing to the second order in the field 



1 



\^ e (ip e )\^ h (ip h )\^ 
J — 2 (rrih cos iph — m e cos p e ) ~^ c n w i p 



4tt 2 

C at XV <j 



(D5) 



+2 (to 2 , cos 2</j h + m 2 e cos 2<p e ) lAk^lp 2 



e 2 R e 



4flM ex COSl^e COS if h - 



CN w j 



F 



where 



tf(i-2) r ^(j-3) 



(.7 -I) 2 - J 2 



(j-2) 2 -j 2 (j + l) 2 -J 2 



[(j + l) 2 -J 2 ][(j + 2) 2 -j 2 ] ■ 

Plugging Eqs. (|D5|) and (001) into Eq. flU]) and noticing 
that the integrals involving linear combinations of the 
cosine-functions are strongly suppressed due to the inte- 
gration over the cosine period, and are therefore negligi- 
ble compared to the one involving the quadratic cosine- 
combination, we obtain 



p2ir p2-rr 
X / dlf e / dlfih- 

Jo JO 



1 — 2 COS (fi e COS l^ftAj (F) 

{z 2 + 4 J R 2 7Jv sin 2 [(^ e -^ l )/2]} 1/2 



(D6) 



with Aj(F) given by Eq. 

The next step is to perform the double integration in 
Eq. (|D6|) . We have to evaluate the two double integrals. 
They are 



h 
and 
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dip, 



2tt 



difh 



p2l! 

h= dip, 
Jo 



o{z 2 + 4i? 2 7Jv sin 2 [(^ e --^)/2]} 1/2 

difh COS ip e COS iph 

o{z 2 + 4i? 2 7W sin 2 [(^ e -^)/2]}V2 



(D7) 



(D8) 



We first notice that both Ii and I2 can be equivalently 
rewritten as follows 



2rr />2tt />2tt rife 

dip A dip h ... = 2 dip e dip h . 
Jo Jo Jo 



(D9) 



due to the symmetry of the integrands with respect to 
the (ip e = iph)-li'ac. Using this property, we substitute iph 
with the new variable t — sin[(p e — iph)/2] in Eqs. (1D7|) 
and (|D8|) . This, after simplifications, yields 



r2ir />sin((/? e /2) 

4 / dip 







dt 



and 



o [(l-t 2 )(z 2 + 4 J R 2 7Ar t 2 )]V2 
dt{l~2t 2 ) 



(D1Q) 



/2 ^ 4 r^ cosV v [(i-^)^ 2 ^^ 2 )]^- 

(Dll) 

Here, the inner integrals are reduced to the incomplete 
elliptical integrals of the first and second kinds (see, e.g., 
Ref. ||). 

We continue the evaluation of Eqs. (|D10[1 and (|D1 1|) 
by expanding the denominators of the integrands in series 
at large and small \z\ as compared to the CN diameter 
2Rcn- One has 



1 



1 



(z 2 + 4^2)1/2 \z\ 



1 ( 2R CN t \- 

2 V M J 



3 { 2R CN t \ 4 _5_ ( 2R CN t \ 6 
8 W 16 I W J 



+ 



for |z|/2i? CJ v > 1, and 



o [(l-< 2 )(z 2 +4i? 2 7 ^ 2 )]V2 2i? CJV 



lim 



»n( V ./S0 /(t) 

at- 



(\z\/2R CN )^0 J\ Z \/2R CN Wl ~ t 2 

for |z|/2i?c , Af "C 1 [f(t) is a polynomial function]. Using 
these in Eqs. ([D10|> and ([DTTj) . we arrive at 



4tt 



4?^ 



i -i(^J+&(^) 4 



2Rcn 



< 1 



2Rcn 



> 1 
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and 



47T 

Rcn 



1 ] J^Rcn \ 



_l+3 ( N_V 
1 + H ^fe ) 



7r 2 /" 2_Rcjv N \ 



, 3 ( 2R0A 
~ 3 \W~) 



' 2i? CJV 

Jfl 
' 2R CN 



<C 1 



> 1 



Plugging these I\ and I2 into Eq. 
leading expansion terms yields 



D6[) and retaining only 



J v 

ireRcN 

M_ 



e 2 » 1 



(D12) 



We see from Eq. (|D12j) that, to the leading order in 
the series expansion parameter, the perpendicular elec- 
trostatic field does not affect the longitudinal electron- 
hole Coulomb potential at large distances \z\ ^> 2Rcn, 
as one would expect. At short distances \z\ -C 2R C n the 
situation is different, however. The potential decreases 
logarithmically with the field dependent amplitude as \z\ 
goes down. The amplitude of the potential decreases 
quadratically as the field increases [see Eq. (|3"3")) ] , thereby 
slowing down the potential fall-off with decreasing \z\, 
or, in other words, making the potential shallower as the 
field increases. Such a behavior can be uniformly approx- 
imated for all \z\ by an appropriately chosen attractive 
cusp-type cutoff potential with the field dependent cut- 
off parameter. Indeed, consider the dimcnsionless func- 
tion f(y) = —2Rcn£ ^eff/e 2 of the dimensionless variable 
y = \z\/2Rcn- Then, according to Eq. (|E)12|) . one has 



$ 1 ( 2 /) = |[l-A J (f)]ln(|), 0<y«l 



$a(l/) = J, y»i 

Now introduce the function 



1 



y + ya 



(D13) 



with the cutoff parameter yo selected in such a way as 
to satisfy the condition $(1) = [$i(l) + $ 2 (l)]/2. This 
yields 



yo = 



7r-21n2[l- Aj(F)} 
•21n2[l- Aj(F)] 



(D14) 



Figure [8] shows the zero-field behavior of the Q(y) 
function as compared to the corresponding < &i(y) and 
4> 2 (y) functions. We see that <&(?/) gradually approaches 
( J ) 2(y) = 1/y f° r increasing y > 1. For decreasing y < 1, 
on the other hand, is very close to the logarith- 

mic behavior as given by <&i(y), with the exception that 
there is no divergence at y ~ due to the presence of the 



— l!y 

(2/x)ln(2/j) 




FIG. 8: (Color online) The dimensionless function (|D13|I with 
the zero-field cutoff parameter (|D14[) . See text for details. 



cutoff. The cutoff parameter (|D14[) is field dependent, 
decreasing as the field grows, which is consistent with 
the behavior of the original potential (|D12j) . Multiply- 
ing Eq. (|D13|) by the dimensional factor —e 2 /2Rcn£ and 
putting y — \z\/2Rcn, we obtain the attractive longitu- 
dinal cusp- type cutoff potential (|31|) we build our analysis 
on in this paper. 



[1] R.Saito, G.Dresselhaus, and M.S.Dresselhaus, Science of 
Fullerens and Carbon Nanotubes (Imperial College Press, 
London, 1998). 

[2] H.Dai, Surf. Sci. 500, 218 (2002). 

[3] L.X.Zheng, M.J.O'Connell, S.K.Doorn, X.Z.Liao, Y.H. 
Zhao, E.A.Akhadov, M.A.Hoffbauer, B.J.Roop, Q.X.Jia, 
R.C.Dye, D.E.Peterson, S.M.Huang, J.Liu and Y.T.Zhu, 
Nature Materials 3, 673 (2004). 

[4] S.M.Huang, B.Maynor, X.Y.Cai, and J.Liu, Advanced 
Materials 15, 1651 (2003). 

[5] R.H.Baughman, A.A.Zakhidov, and W.A.de Heer, Sci- 
ence 297, 787 (2002). 



[6] A.Popescu, L.M.Woods, and I.V.Bondarev, Nanotech- 

nology 19, 435702 (2008). 
[7] J.E.Trancik, S.C.Barton, and J.Hone, Nano Lett. 8, 982 

(2008). 

[8] I.V.Bondarev and B.Vlahovic, Phys. Rev. B 74, 073401 

(2006) . 

[9] I.V.Bondarev and B.Vlahovic, Phys. Rev. B 75, 033402 

(2007) . 

[10] I.V.Bondarev, J. Electron. Mater. 36, 1579(2007). 
[11] I.V.Bondarev, Optics & Spectroscopy (Springer) 103, 
366 (2007). 

[12] H.Shimoda, B.Gao, X.P.Tang, A.Kleinhammes, L.Flc- 



21 



ming, Y.Wu, and O.Zhou, Phys. Rev. Lett. 88, 015502 
(2001). [43 
G.-H. Jeong, A.A.Farajian, R.Hatakeyama, T.Hirata, 
T.Yaguchi, K.Tohji, H.Mizuseki, and Y.Kawazoe, Phys. [44 
Rev. B 68, 075410 (2003) 

G. -H. Jeong, A.A.Farajian, T.Hirata, R.Hatakeyama, K. [45 
Tohji, T.M.Briere, H.Mizuseki, and Y.Kawazoe, Thin 

Solid Films 435, 307 (2003). [46 
M.Khazaei, A.A.Farajian, G.-H. Jeong, H.Mizuseki, T. 
Hirata, R.Hatakeyama, and Y.Kawazoe, J. Phys. Chem. [47 
B 108, 15529 (2004). 

A.Hogele, Ch.Galland, M. Winger, and A.Imamoglu, 

Phys. Rev. Lett. 100, 217401 (2008). [48 

T.G.Pedersen, Phys. Rev. B 67, 073401 (2003). 

T.G.Pedersen, Carbon 42, 1007 (2004). [49 

R.B.Capaz, C.D.Spataru, S.Ismail-Beigi, and S.G.Louie, 

Phys. Rev. B 74, 121401(R) (2006). [50 

C.D.Spataru, S.Ismail-Beigi, L.X.Benedict, and S.G. 

Louie, Phys. Rev. Lett. 92, 077402 (2004). [51 

F.Wang, G.Dukovic, L.E.Brus, and T.F. Heinz, Phys. 

Rev. Lett. 92, 177401 (2004). 

F.Wang, G.Dukovic, L.E.Brus, and T.F. Heinz, Science [52 
308, 838 (2005). 

A. Hagen, M.Steiner, M.B.Raschke, C.Lienau, T.Hertel, [53 

H. Qian, A.J.Meixner, and A.Hartschuh, Phys. Rev. Lett. [54 
95, 197401 (2005). 
F.Plentz, H.B.Ribeiro, A.Jorio, M.S.Strano, M.A.Pimen- [55 
ta, Phys. Rev. Lett. 95, 247401 (2005). 

F.Xia, M.Steiner, Y. -M.Lin, and Ph.Avouris, Nature [56 
Nanotechnology 6, 609 (2008). [57 

B. F.Habenicht and O.V.Prezhdo, Phys. Rev. Lett. 100, 
197402 (2008). [58 
V.Perebeinos, J.Tersoff, and Ph.Avouris, Phys. Rev. 

Lett. 94, 027402 (2005). 

M.Lazzeri, S.Piscanec, F.Mauri, A.C.Ferrari, and J.Ro- [59 
bertson, Phys. Rev. Lett. 95, 236802 (2005). 
S.Piscanec, M.Lazzeri, J.Robertson, A.C.Ferrari, and F. [60 
Mauri, Phys. Rev. B 75, 035427 (2007). 

S.Zaric, G.N.Ostojic, J. Shaver, J.Kono, O.Portugall, [61 
P.H.Frings, G.L.J.A.Rikken, M.Furis, S.A.Crooker, X. [62 
Wei, V.C.Moore, R.H.Hauge, and R.E.Smalley, Phys. 
Rev. Lett. 96, 016406 (2006). [63 
V.Perebeinos and Ph.Avouris, Nano Lett. 7, 609 (2007). 
A.Srivastava, H.Htoon, V.I.Klimov, and J.Kono, Phys. [64 
Rev. Lett. 101, 087402 (2008). 

M.S.Dresselhaus, G.Dresselhaus, R.Saito, and A.Jorio, [65 
Annu. Rev. Phys. Chem. 58, 719 (2007). [66 
T.Pichler, M.Knupfer, M.S. Golden, J.Fink, A.Rinzler, [67 
and R.E.Smalley, Phys. Rev. Lett. 80, 4729 (1998). 

C. D.Spataru, S.Ismail-Beigi, R.B.Capaz, and S.G.Louie, 
Phys. Rev. Lett. 95, 247402 (2005). [68 
Y.-Z.Ma, C.D.Spataru, L.Valkunas, S.G.Louie, and G.R. 
Fleming, Phys. Rev. B 74, 085402 (2006). 

J.Bellessa, C.Bonnand, J.C.Plenet, and J.Mugnier, Phys. [69 
Rev. Lett. 93, 036404 (2004). 

W.Zhang, A.O.Govorov, and G.W.Bryant, Phys. Rev. 
Lett. 97, 146804 (2006). [70 
Y.Fedutik, V.V.Temnov, O.Schops, U.Woggon, and M.V. 
Artemyev, Phys. Rev. Lett. 99, 136802 (2007). 

I. V.Bondarev, G.Ya.Slepyan and S.A.Maksimenko, Phys. [71 
Rev. Lett. 89, 115504 (2002). 

I.V.Bondarev and Ph.Lambin, Phys. Rev. B 70, 035407 
(2004). [72 
I.V.Bondarev and Ph.Lambin, Phys. Lett. A 328, 235 



(2004) . 

I.V.Bondarev and Ph.Lambin, Solid State Commun. 
132, 203 (2004). 

I.V.Bondarev and Ph.Lambin, Phys. Rev. B 72, 035451 

(2005) . 

I.V.Bondarev and Ph.Lambin, in: Trends in Nanotubes 
Reasearch (Nova Science, NY, 2006). Ch.6, pp. 139-183. 
W.Vogel and D.-G.Welsch, Quantum Optics (Wiley- 
VCH, NY, 2006). Ch.10, p.337. 

L. Knoll, S.Scheel, and D.-G.Welsch, in: Coherence and 
Statistics of Photons and Atoms, edited by J.Pefina (Wi- 
ley, NY, 2001). 

S.Y.Buhmann and D.-G.Welsch, Prog. Quantum Elec- 
tron. 31, 51 (2007). 

L.X.Benedict, S.G.Louie, and M.L.Cohen, Phys. Rev. B 
52, 8541 (1995). 

S.Tasaki, K.Maekawa, and TYamabe, Phys. Rev. B 57, 
9301 (1998). 

Z.M.Li, Z.K.Tang, H.J.Liu, N.Wang, C.T.Chan, R.Saito, 
S.Okada, G.D.Li, J.S.Chen, N.Nagasawa, and S.Tsuda, 
Phys. Rev. Lett. 87, 127401 (2001). 

A. G.Marinopoulos, L. Reining, A.Rubio, and N.Vast, 
Phys. Rev. Lett. 91, 046402 (2003). 

TAndo, J. Phys. Soc. Jpn. 74, 777 (2005). 

B. Kozinsky and N.Marzari, Phys. Rev. Lett. 96, 166801 

(2006) . 

H.Haken, Quantum Field Theory of Solids, (North- 
Holland, Amsterdam, 1976). 

S.Uryu and TAndo, Phys. Rev. B 74, 155411 (2006). 
L.D. Landau and E.M.Lifshits, The Classical Theory of 
Fields (Pergamon, New York, 1975). 
K.Kempa and P.R.Chura, special edition of the Kluwer 
Academic Press Journal, edited by L.Liz-Marzan and 
M.Giersig (2002). 

K.Kempa, D.A.Broido, C.Beckwith, and J.Cen, Phys. 
Rev. B 40, 8385 (1989). 

A.S.Davydov, Quantum Mechanics (Pergamon, New 
York, 1976). 

V.N.Popov, L.Henrard, Phys. Rev. B 70, 115407 (2004). 
M.F.Lin, D.S.Chuu, and K.W.-K.Shung, Phys. Rev. B 
56, 1430 (1997). 

H.Ehrenreich and M.H.Cohen, Phys. Rev. 115, 786 
(1959). 

X. Blase, L.X.Benedict, E.L.Shirley, and S.G.Louie, Phys. 
Rev. Lett. 72, 1878 (1994). 
K.Kempa, Phys. Rev. B 66, 195406 (2002). 
E.Hanamura, Phys. Rev. B 38, 1228 (1988). 
J.P.Reithmaier, G.S§k, A.Lomer, C.Hofmann, S.Kuhn, 
S.Reitzenstein, L.V.Keldysh, V.D.Kulakovskii, T.L.Rei- 
necke, and A.Forchel, Nature 432, 197 (2004). 
TYoshie, A.Scherer, J.Hendrickson, G.Khitrova, H.M. 
Gibbs, G.Rupper, C.E11, O.B.Shchekin, and D.G.Deppe, 
Nature 432, 200 (2004). 

E.Peter, P.Senellart, D.Martrou, A.Lemaitre, J. Hours, 
J.M.Gerard, and J.Bloch, Phys. Rev. Lett. 95, 067401 
(2005). 

D.A. B.Miller, D.S.Chemla, T.C.Damen, A.C.Gossard, 
WWiegmann, T.H.Wood, and C.A.Burrus, Phys. Rev. 
Lett. 53, 2173 (1984). 

D.A. B.Miller, D.S.Chemla, T.C.Damen, A.C.Gossard, 
W.Wiegmann, T.H.Wood, and C.A.Burrus, Phys. Rev. 
B 32, 1043 (1985). 

A.Zrenner, E.Beham, S.Stufler, F.Findeis, M.Bichlcr, 
and G.Abstreiter, Nature 418, 612 (2002). 



22 



[73] H.J.Krenner, E.C.Clark, T.Nakaoka, M.Bichler, C.Sche- 
urer, G.Abstreiter, and J.J.Finley, Phys. Rev. Lett. 97, 
076403 (2006). 

[74] In real CNs, the existence of two equivalent energy valleys 
in the 1st Brillouin zone, the K- and Tf'-valleys with 
opposite electron helicities about the CN axis, results 
into dark and bright excitonic states in the lowest energy 
spin-singlet manifold [75[. Since the electric interaction 
does not involve spin variables, both K- and if'-valleys 
are affected equally by the electrostatic field in our case, 
and the detailed structure of the exciton wave function 
multiplet is not important. This is opposite to the non- 
zero magnetostatic field case where the field affects the 
K- and if '-valleys differently either to brighten the dark 
excitonic states [3^| , or to create Landau sublevels [53| for 
longitudinal and perpendicular orientation, respectively. 

[75] T.Ando, J. Phys. Soc. Jpn. 75, 024707 (2006). 

[76] T.Ogawa and T.Takagahara, Phys. Rev. B 44, 8138 
(1991). 

[77] A.Jorio, C.Fantini, M.A.Pimenta, R.B.Capaz, Ge.G.Sa- 
msonidze, G.Dresselhaus, M.S.Dresselhaus, J.Jiang, N. 
Kobayashi, A.Griineis, and R.Saito, Phys. Rev. B 71, 
075401 (2005). 

[78] We are based on the zero-exciton-temperature approxi- 
mation in here [7j|, which is well justified because of the 
exciton excitation energies much larger than ksT in CNs. 
The exciton Hamiltonian Q does not require the thermal 
averaging over the exciton degrees of freedom then, yield- 
ing the temperature independent total exciton energy (0 . 
One has to keep in mind, however, that the exciton exci- 



tation energy can be affected by the enviromental effect 

not under consideration in here (see Ref. [80l]V 
[79] A.Suna, Phys. Rev. 135, Alll (1964). 
[80] Y.Miyauchi, R.Saito, K.Sato, Y.Ohno, R.Iwasaki, T.Mi- 

zutani, J.Jiang, and S.Maruyama, Chem. Phys. Lett. 

442, 394 (2007). 
[81] I.V.Bondarev, Y.Nagai, M.Kakimoto, T.Hyodo, Phys. 

Rev. B 72, 012303 (2005). 
[82] J.D.Jackson, Classical Electrodynamics (Wiley, New 

York, 1975). 

[83] Higher Transcendental Functions, edited by H.Bateman 

and A.Erdelyi (Mc Graw-Hill, New York, 1955). Vol. 3. 
[84] Handbook of Mathematical Functions, edited by M.Ab- 

ramovitz and I.A.Stegun (Dover, New York, 1972). 
[85] L.D. Landau and E.M.Lifshits, Quantum Mechanics: 

Non-Relativistic Theory (Pergamon, New York, 1977). 
[86] S.Liu, J.Li, Q.Shen, Y.Cao, X.Guo, G.Zhang, C.Feng, 

J.Zhang, Z.Liu, M.L.Steigerwald, D.Xu, and C.Nuckolls, 

Angew. Chem. 48, 4759 (2009). 
[87] P.L.Hernandes-Martinez and A.O.Govorov, Phys. Rev. B 

78, 035314 (2008). 
[88] H.Qian, C.Georgi, N.Anderson, A.A.Green, M.C.Her- 

sam, L.Novotny, and A.Hartschuh, Nano Lett. 8, 1363 

(2008). 

[89] P.H.Tan, A.G.Rozhin, T.Hasan, P.Hu, V.Scardaci, 
W.I.Milne, and A.C.Ferrari, Phys. Rev. Lett. 99, 137402 
(2007). 

[90] H.Hirori, K.Matsuda, and Y.Kanemitsu, Phys. Rev. B 
78, 113409 (2008). 



